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We investigate the influence of inhomogeneity in the pairing coupling constant U(r) on dirty BCS 
superconductors, focusing on T c , the order parameter A(r), and the energy gap E g (r). Within 
mean-field theory, we find that when the length-scale of the inhomogeneity is comparable to, or 
larger than the coherence length, the ratio 2E g /T c is significantly reduced from that of a homoge- 
neous superconductor, while in the opposite limit this ratio stays unmodified. In two dimensions, 
when strong phase fluctuations are included, the Kosterlitz-Thouless temperature Tkt is also stud- 
ied. We find that when the inhomogeneity length scale is much larger than the coherence length, 
2E g /TpcT can be larger than the usual BCS value. We use our results to qualitatively explain recent 
experimental observation of a surprisingly low value of 2E g /T c in thin films. 



I. INTRODUCTION 

The presence of disorder in essentially all supercon- 
ducting systems makes research of the interplay of disor- 
der and superconductivity essential. In their pioneering 
work, Anderson^, and Abrikosov and Gorko\*2, claimed 
that nonmagnetic impurities have no considerable effect 
on the thermodynamic properties of s-wave supercon- 
ductors; this result is known as "Anderson theorem" for 
weakly disordered dirty superconductors. Since the dis- 
covery and elucidation of the localization phenomenon 3 , 
corrections to the Anderson theorem have been in- 
tensively investigated both experimentally^'^' 7 -' 8 - and 
theoretically^Oaiii^^iii^iieaiiis. Within mean field 
theory, it has been shown that if one neglects Coulomb in- 
teractions, pairing survives below the mobility edge until 
the localization length reaches a critical valu o 11 ' 12 . But 
interactions change this picture significantly, since the ef- 
fect of Coulomb repulsion is strengthened by localization, 
resulting in a suppressed effective attractive interaction 
and thus a reduced mean-field r^ i 10 i 13 i 14 i 15 . An under- 
lying assumption of these works is the uniformity of the 
superconducting order parameter, which has been ques- 
tioned by numerical simulations in recent year a 17 ' 18 . 

Experiments in this field focused on two-dimensional 
(2d) superconductors, namely superconducting thin 
films. The disorder in superconducting films is expected 
to reduce the superfluid density and the phase order- 
ing temperature, i.e., the Kosterlitz-Thouless temper- 
ature Tkt, in addition to suppressing the mean field 
T c . These considerations naturally lead to the pos- 
sibility of a quantum superconductor-insulator transi- 
tion (SIT) at a critical amount of disorder or mag- 
netic field. Furthermore, the scale invariant na- 
ture of a film's resistance raised expectations that 
such an SIT would exhibit many universal features^ 9 -. 
The superconducting-insulator transition was intensively 
studied experimentally 20 . 21 i 22 i 23 i 24 . 25 i 26 i 27 i 28 i 29 i 30 . 31 . The 
theoretical viewpoint on these transitions took two main 
forms: the nature of the SIT was interpreted either 
as the breaking of Cooper pairs caused by amplitude 
fluctuatio n 14 i 15 i 16 , or localization of Cooper pairs re- 



sulting from phase fluctuationi^^i^ 3 .^^! 3 ^. While the 
nature of the SIT in various systems is still debated, 
in recent years the interest in this problem is further 
intensified by the observation of a possible metallic 
phase intervening the superconducting and insulating 
phas o 37 i 38 ' 39 i 40 ' 41 . This observation stimulated several 
theoretical proposal o 42 i 43 ' 44 i 45 ' 46 i 47 , but its origin is still 
a mystery. 

Motivated by the thin-film physics, more experimen- 
tal studies focused on the nature of the density of states 
(DOS) and the quasi-particle energy gap of disordered 
single layer superconducting thin films^ ' 7 ' 23 ' 25 ' 26 ' 28 and 
superconductor - normal-metal (SN) bilayer a 41 ' 48 ! 49 . In- 
terestingly, these studies found a broadening of the BCS 
peak and also a subgap density of states^' 2 ^'^' 2 ^'^ 8 .^ 8 -. 
Of particular interest to us is the work in Ref. H^, which 
studied a thin SN bilayer system, and found a surpris- 
ingly low value of the ratio of the energy gap to T c , in 
contradiction to standard BCS theory, and the theory of 



proximity 



50,51.52 



where it is claimed that the energy gap- 



T c ratio should be bounded from below by ~ 3.52. A drop 
below this bound, 2E g /T c < 3.52, was also observed in 
amorphous Bi films as it approaches the disorder tuned 
SI T 23 ' 26 . Similar trends were also observed in SN bilayers 
in Ref. Hi] and in amorphous tin films in Ref. 0. 

In this paper we show that a reduction of the 2E g /T c 
ratio in a dirty superconductor could be explained as a 
consequence of inhomogeneity in the pairing interaction. 
In SN bilayer thin films, thickness fluctuations of either 
layer result in effective pairing inhomogeneity (in thin 
SN bilayers the effective pairing is the volume averaged 
one, c.f., Ref. I5l"lf52l and Sec. llVj) . Such inhomogeneities 
in other systems occur due to grain boundaries, dislo- 
cations, or compositional heterogeneity in alloys^ 3 -. For 
simplicity we will assume in our analysis that the pair- 
ing coupling constant takes a one-dimensional modulat- 
ing form: 



U{r) = U + U Q cos(Qx). 



(1) 



In bilayer SN films, the effect of localization and 
Coulomb interaction is minor compared to proximity ef- 
fect, and therefore we will neglect these complications in 
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this work. 



In our results, the ratio between the inhomogeneity 
length, L = l/Q, and the superconducting coherence 
length £, plays a crucial role. When Q£ ^> 1, the su- 
perconducting properties are determined by an effective 
coupling U < U e ff < U + Uq 54 . In this limit, the ratio 
2E g /T c is preserved at the standard BCS value ~ 3.52. 
Small corrections are obtained when 1/(Q£) is finite. In 
the opposite limit, Q( < 1, the system tends to be de- 
termined by the local value of U(x). Within mean field 
theory, the ratio 2E g /T c is generally suppressed from 
the BCS value 3.52; in 2d, however, when one includes 
the thermal phase fluctuation and studies the Kosterlitz- 
Thouless temperature, Tkt-, the ratio 2E 9 /Tkt can be 
larger than the usual BCS value. These results on 2E g /T c 
are summarized in FIG. [6] 



Our analysis is inspired by similar previously stud- 
ied models. Particularly, the T c of the clean case of 
this model has been analyzed in Ref. 54. Here we 
extend the study of non-uniform pairing to both T c 
and zero-temperature properties of disordered films, in 
the regime where the electron mean free path I obeys 
1/kp <C I <C £o ~ which is relevant to the ex- 

periments of Long et al . 48 ' 49 . Note that while Anderson 
theorem states that the critical temperature and gap of a 
homogenous superconductor do not depend on disorder—, 
in an inhomogeneous system the theorem does not hold. 
Indeed, we find that the results of Ref. are modi- 
fied in the dirty case. In another related work, a system 
with a Gaussian distribution of the inverse pairing inter- 
action was studie d 55 ' 56 . It was shown that an exponen- 
tially decaying subgap density of states appears due to 
mesoscopic fluctuations which lie beyond the mean field 
picture. Finally, inhomogeneous coupling in the attrac- 
tive Hubbard mode& and lattice XY model^ 8 - were also 
analyzed, with relevance to High-T c materials. 



This paper is organized as follows. In Sec. |TT] we re- 
view the quasiclassical Green's function formalism which 
we use, and briefly demonstrate how it works for the 
usual dirty superconductors with spatially uniform cou- 
pling constant. Then, in Sec. IIIII we discuss the cases 
with nonuniform coupling classified by the competition of 
two length scales: the coherence length £ and the length 
scale associated with the variation of the coupling con- 
stant L — l/Q. We will also discuss the effect of other 
types of inhomogeneities briefly. In section IIVI we pro- 
vide a useful analogy with superconductor-normal metal 
superlattice to provide more physical intuition about our 
results on the energy gaps. In section [V] we will summa- 
rize our analysis and discuss the connection with experi- 
mental results. 



II. THE GAP EQUATION OF A NONUNIFORM 
FILM 

The starting point of our analysis is the standard s- 
wave BCS Hamiltonian: 

H = Hq + Hi n t + Himp, 

Hint = -C/(r)^(r)^(r)^ T (f)^(r), (2) 

where £ = — ^ — A 4 , and U{r) > is the attractive 
coupling constant between electrons, and Hi mp includes 
scattering with nonmagnetic impurities. When the pair- 
ing interaction, U(r), is nonuniform, so is the order pa- 
rameter in this system. A standard technique to tackle 
this non-uniform superconductivity problem is the qua- 
siclassical Green's function o 59 ' 60 ' 61 . In the dirty limit 
t <C £o ~ y-j the quasiclassical Green's functions obey a 
simple form of the Usadel equation, which in the absence 
of a phase gradient is: 

— (-V 2 (9) = Acos(9-w„sin6>, (3) 

where D = 2 v fI is the diffusion constant, I is the mean 
free path, d is the spatial dimension, and A is the super- 
conducting order parameter. 9 is a real function of space 
and Matsubara frequencies u n and is a parametrization 
of the quasiclassical Green functions g and /: 

g = cos 9, f = p = -i sin 9. (4) 

Also, we list the relation between the integrated quasi- 
classical Green's function and Gor'kov's Green's function 
G and F: 

«>- f I 

where r is the center of mass coordinate, and p is mo- 
mentum corresponding to the relative coordinate; £l p is 
the angle of momentum p and Np is the density of states 
(per spin) of the normal state at the Fermi energy. The 
self-consistency equation reads: 

A(r) = U(f)N F 7rTj2iUAr)- (5) 

n 

For simplicity we assume the pairing is as given in Eq. 

©, 

U(r) = U + U Q cos(Qx). 
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A. The uniform pairing case 

Before analyzing the inhomogeneous pairing problem, 
let us briefly review the calculation of T c , the supercon- 
ducting order parameter A(T = 0), and the DOS v{E) 
of a dirty superconductor with a spatially uniform cou- 
pling constant U , using quasiclassical Green's functions. 
In this case Eqs. (j3|) and ([5]) admit a uniform solution 
for both 9 and A: 



9 = arctan 



A 



(6) 



Using ©, we obtain the standard BCS self-consistency 
equation: 



1 = UN F nT V —=L 



T c and A(T = 0) are easily obtained from (JT]): 



(7) 



2C 



-uj D e 



i \t=o) 



where C = e 7 « 1.78, with 7 = 0.5772 . . . the Euler con- 
stant, and u>d the Debye frequency. The DOS can be 
obtained from the retarded quasiclassical Green's func- 
tion: u(E) = Re{g R (E)}, which can be obtained from 
g(uj n ) = cos(0 n ) by analytical continuation iu> — > E+i0 + : 



v(E) = Re. 



-%E 



^/A 2 - {E + i0+) 2 



0, if E < A 



Thus there exists a gap in the excitation spectrum E g — 
A, and its ratio with T c is a universal number 7r/C ~ 
1.76. As expected, these results for dirty superconductors 
are exactly the same as those of clean superconductors, 
thus explicitly illustrating Anderson theorem. 



III. THE CASE OF INHOMOGENEOUS 
PAIRING 

Using the formalism reviewed in the previous section, 
we now discuss the non-uniform superconducting film. 
Our discussion will concentrate on the limits of fast and 
slow pairing modulations, i.e., large and small Q£ respec- 
tively (£ is the zero temperature coherence length in the 
dirty limit: £ = y/hD /A T =o ~ y/hD/T c , where A is the 
spatially averaged A(x)). 



A. Fast pairing modulation: proximity enhanced 
superconductivity 

With a nonuniform coupling U(x), uniform solution of 
either 9{x) or A(x) no longer exists. When fast pairing 
modulation are present, the angle 9 is dominated by its 
k = Fourier component, 9q, since it can not respond 



faster than its characteristic length scale £. Corrections 
to the uniform solution are of the form 9\ cos(Qx), and 
are suppressed by powers of From Eq. {I}, we see 
that in contrast to 9, the order parameter A(x) has a 
factor of U(x) in its definition, and therefore it can fluc- 
tuate with the fast modulation of U (x) . The modulating 
component of A(x) is thus only suppressed by Uq/U, 
while the modulating part of 6{x) is suppressed by both 
Uq/U and 1/(Q£). Keeping both 1/Q£ <C 1 and expand- 
ing in Uq/U, we can perturbatively solve Eqs. © and 
fl5j). Starting with: 

A(x) = A + Ai cos(Qx),9(x) = 9 + 9x cos(Qx); (8) 

Eq. ([3|) can be solved order by order: 

'A 



arctan 



Ai 



(9) 



fQ^V^+Ag + ^ + Ag 



The self-consistency equation ([5]) can be Fourier trans- 
formed: 



A = N F nT [ U sin 9 + 2 



, Uq cos go 
2 2 



, (10) 



£ = N F ,Tj2[u^9 1 + ^ S ine ), 



cos9 U Q . 



where the u) n index of 9q and 9\ is implicit. 

When T — > T c , we can linearize #0 and #1 with respect 
to Aq and Ai, respectively: 



sm9 ~-^-r, 6»i(cos6» ) ~ — ^ 



-DQ 2 
2 



Note that 
Wo 1 

V » In JV + 21n2 + 7 for 7V » 1, (11) 

z — ' n + 1/2 



n=0 



where 7 is the Euler constant, we have approximately 

2ttT V — « ln^CWvrT), (12) 

2ttT 



w„=0 

1 



« n + £Q 2 /2 



In 1 



DQZ/2) ' 



where, as before, C = e 7 w 1.78 and wd is the Debye 
frequency. Defining 



JT = £W F ln(2Co; D /7rT), #1 = IW F i n 1 



2w z 
£Q 2 



(13) 
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we get 

Ai = ^S-KqAq + KxAx. 

T c is the temperature at which this equation admits a 
nonzero solution: 



2C / 1 
T c = — exp - — . , 

7T V U eff N F J 



where the effective pairing strength is: 



V.ss ' \ i l'^ 



[7 J 2(1- Kx) 



(14) 



(15) 



This is the dirty case analogue of the result obtained by 
Ref. 

Next we turn to the order parameter. At T = the 
sums in the self-consistency equations pU|) become inte- 
grals, which can be performed (see also Appendix [A"| : 



A = N F U A \n 

Ai = KxAx 
2 2 

thus giving the solution 



2uJl 
A 
N F U Q 



A 0(T=0) = 2uj d exp I - 



-"(If) 



1 



(16) 



Ai 



(T=0) 



*0(T=0) 



UeffN F J' 

Uq 1 



C7, 



with the same U e ff defined in (fl5j) . Noting that Ao is 
the spatially averaged value of the order parameter A, 
we arrive at the conclusion that in the limit Qt; S> 1, the 
ratio 



2A 

~t7 



2A 



2n 
~C 



(17) 



is preserved. 

The modification of the gap, however, must be ad- 
dressed separately. Although the gap and the order pa- 
rameter coincide for a uniform BCS superconductor, this 
is not generally true in a nonuniform superconductor. 
To obtain the DOS and the gap one has to rephrase the 
problem in a real-time formalism and calculate the re- 
tarded Green's function which is parameterized by a com- 
plex 9(x, E) = 6'(x, E) + i9"(x, E) with both 0', 9" real, 
and then compute the DOS via v(x, E) = Keg R (x, E) = 
Re cos 9 (x, E) = cos 6' cosh g / /60 i 61 . Naively one can per- 
form the prescription iuj — > E+iO + in the imaginary time 
Green's functions to obtain the retarded ones, but our 
perturbative solution will break down as E approaches 
Ao, since 9\ diverges faster than 9q. Therefore to ana- 
lyze the gap one has to re-solve the real time counterpart 
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FIG. 1: 

Q£ for Q£ 3> 1. The two curves are for Ai/Ao = 0.1 and 
0.2, respectively. Here, Ao and Ai are the uniform and os- 
cillating components of the order parameter, respectively. Q 
is the modulating wavevector of the inhomogeneous coupling 
constant; £ is the superconducting coherence length. The esti- 
mated numerical error of E g j Ao is about 0.01. The deviation 
of Eg from Ao is small, but it increases with larger Ai/Ao or 
smaller Q£. 



of equation ([3]) with A(a;) given above. Note that our so- 
lution of A[x) is still valid, sparing us the need to solve 
the self-consistency equation. 
In real time, Eq. becomes: 

-^d 2 J' = cos0'(Acosh0"-£sinh0"), 



D 



d 2 x 9" = sin0'(Asinh0" - £cosh0"). (18) 



We numerically solved these coupled equations with pe- 
riodic boundary condition on [0, 2ir/Q], and computed 
the DOS v(E) = cos0i cosh 02, and thereby obtained the 
gap. We find that despite the fluctuating A(x), the en- 
ergy gap, Eg, is spatially uniform. Fig. [T] shows a graph 
of Eg vs. Qt; for Ai/Ao = 0.1 and 0.2. Again, in the plot 
we define the coherence length £ to be ^/fiD/A^o = 
^JfiD / Ao i t=o- One can see that in the limit Qt; — > oo 
Eg coincides with Ao, and nonzero 1/(Q£) brings about 
only small corrections to make the gap slightly smaller 
than Ao. These corrections increase with smaller Qt; or 
larger U Q /U (i.e., Ai/A ). Thus we find that for Qt; > 1 



IE, 



< 



2A 



(T=0) 



T- 



2tt 

= c = 3 - 52 ' 



(19) 



It is easy to understand the uniformity of E g , since the 
wave function of a quasiparticle excitation should be ex- 
tended on a length scale 1/Q <C £. Some intuition for 



the fact that E g « Aq is provided in Sec. [TV. 
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B. Slow pairing fluctuations: WKB-like local 
superconductivity 

When the pairing strength fluctuates slowly, i.e., over a 
large distance, both the Green's functions and the order 
parameter A(x) can vary on the length scale of 1/Q, and 
we can approximate the zeroth order solution by a 'local 
solution': 



9q{x) = arctan 



A(x) 



(20) 



where A(x) is to be solved from the self-consistency equa- 
tion. This 'local' property of the system implies a large 
spatial variation of both A(x) and 6{x), in contrast to the 
Q£ 1 case. To improve the zeroth order solution, we 
write 0(x) = 9q{x~)+6\(x). Neglecting the small gradient 
term of 6\, one can solve for 9\ from Usadel's equation 



D 



co n d 2 x A 



(A 2 + W 2)3/2 (A2 +W 2)5/2 



(21) 



thus the self-consistency equation |[SJ) becomes 

A cj„ 



A(x) = U(x)N f 2ttT 



2-kT 



n=0 



7a 5 



In the Ginzburg-Landau regime, one is justified in keep- 
ing lowest order terms in (|22[) : 



A(x) = U(x)N F \ A(x) In 



'2Cu L 
V ttT 



7C(3) 



8vr 2 T 2 



8T 



A 3 (x) 
9, 2 A(x)l(23) 



where C( n ) is the Riemann £ function. Remarkably, equa- 
tion (|23[) is nothing but the Ginzburg-Landau equation 
for a modulating coupling constant U(x) with Qt; <C 1, 
and is pre cisely the dirty case analogue of equation (9) 
in Ref. 1 54. with £ replaced by the dirty limit expression 
£ 2 = hTrD/8T (£ is slightly different from the coherence 
length defined in this work £ = y/hD/A.T=0i where A 
is the spatially averaged A(x)). In the limit Q£ — * 0, 
A(x) would be determined only by the local value of 
U(x), and the mean field transition temperature would 
be given by T c ^ max = 2Cu> D /Trexp(-l/(U + Uq)). A 
small but nonzero Q£ leads to a weak coupling between 
spatial regions, hence slightly reducing the mean field T c . 
Following the analysis of Ref. one obtains the mean 
field transition temperature: 



T, 



MF 



2Cuj 



D e -l/N F (U+U Q ) e -iQA/V2 



(24) 



where A = ^U Q /{N F U 2 ). 

Although the inhomogeneous U(x) largely increases 
the mean field T c , it also makes the system more sus- 
ceptible to phase fluctuations. This effect will be more 
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FIG. 2: The mean field transition temperature T c , the 
Kosterlitz-Thouless temperature Tkt, and the minimum 
mean field transition temperature T c ,min (see below) (a) vs. 
Uq/U with QJ = 0.3; (b) vs. Uq/U with Q£ = 0.1; (c) vs. 
Qi with Uq/U = 0.1. In all cases UN F = 0.2. T is in units 



of T c . 



1/N F U 



Dashed curve is T c determined 



from equation (|24p : dash-dotted curve is the minimum mean 
field T c given by T c , mln = 2£cj De - 1/NF( -°- ]u <3 l) ; solid curve 
is Tkt obtained from numerically minimizing (125[) and then 
solving (|27|l . Here, U and Uq are the uniform and oscillat- 
ing components of the coupling constant, respectively. N F is 
the density of states of the normal state; Q is the modulating 
wavevector of the inhomogeneous coupling constant; £ is the 
superconducting coherence length. The estimated numerical 
error of Tkt /T c , a is about 0.01. 
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pronounced in a two-dimensional superconductor, which 
we will focus on now. A film becomes superconducting 
through a Kosterlitz-Thouless transition. To determine 
the Kosterlitz-Thouless transition temperature, Tkt, we 
note that the Ginzburg-Landau free energy correspond- 
ing to ([231) is 



F(A(x)) 

a[x) 




AY / d 3 x{a(x)A 2 {x) + ^A 4 (x) 



l(d x A) 2 )} 
1 



(25) 



N F U(x) 
7C(3) 



-In 



2 x 1 



78ujd\ 



nT J 



8tt 2 T 2 ' 



7 



8T ' 



As a functional of A(x), F can be minimized nu- 
merically, thus giving a solution of A(x). The free en- 
ergy cost for phase fluctuations is approximately F = 
\ J d 2 xJ{x)(V0) 2 . For quasi-2d films, 



J(x) = 2N ± Np d £ 2 \A MF (x)\ 



(26) 



where Np is the 2d electron DOS, N±_ is the number of 



TTtlD 



and Amf is the mean field solution 



channels, £ = y -gy - , ana l±mf 
of (|25[) . To explain the bilayer thin film experiments in- 
vestigated by Long et a L 48 ' 49 , we use the measured value 
of the diffusion constant D = 5 x 10 _3 m 2 s~ 1 (see Ref. 
HH), and estimate N± = kpd/n m 50, where the film 
thickness d w 10 ~ 20n m 48 i 49 , and the Fermi wave vec- 
tor kp ~ lA . As in Ref. 54, one can estimate Tkt 
self-consistently from 



T, 



KT 



J(x)(l/J(x))- 



(27) 



since J(x) is the stiffness along the "stripes", while 
(1/ J(x))~ 1 perpendicular to the "stripes" . Although our 
estimation of N± is crude, the value of Tkt is very insen- 
sitive to it. This is because Tkt is solved self-consistently 
from (|2"7| . If one attempts to use a larger N± in (|26)l . 
the enhancement of Tkt is limited by J{x) which itself 
is suppressed as temperature increases. Typical solutions 
of Tkt are shown in FIG. O One can see that the phase 
fluctuation region, i.e. the difference between and 
Tkt, increases with stronger inhomogeneity (FIG. 12a) 
and (b)). Also for longer wave length modulation Tkt is 
reduced more strongly (FIG. Htc)). Heuristically, this is 
because for smaller Q£ the superconducting stripes be- 
come farther apart, and therefore it is more difficult for 
them to achieve phase coherence. 

Moving our focus to the zero-temperature order pa- 
rameter and gap, we note that at T = the integrals in 
equation (j2"2"|) can be done: 



■kD8 2 A vD(d x Af 



AW -A(*)h/ 2 ^ 



U{x)N t 



This can be approximately solved by: 

A(a:) « A (x)e-^\ 
A (x) = 2w D e~ JV A a => , 

T)(x) 



(28) 



A(x) J 8A(x) 16A 2 (x) 



%^Q 2 A 2 (cos(Qx) - l -A 2 s\n 2 {Qx) 



Note that d x A(x) « -A 2 Qsuv(Qx)A(x) [with A de- 
fined under Eq. ([24]) ]. for our WKB analysis to be self- 
consistent, we need to require the that A < 0(1), thus 
Uq/U needs to be small. Also, when this is satisfied, 
rj(x) leads to a slight averaging between A(x), which is 
a manifestation of proximity effect. 

To analyze the gap, we must switch to a real time 
formalism again, since our perturbative solution for the 
Green's function becomes invalid as E — > A(x). Thus 
we have to solve the real time Usadel equation (fP8|) with 
A(x) obtained above. Using the same numerical code as 
in Sec. IIII A[ we have obtained the local gap E g (x), which 
is plotted vs. x in FIG.[3]for half a period of modulation. 
One can see that in general E g (x) is lower than A(x), 
and when Q£ = 0.3, E g (x) is largely set by the region 
with weakest coupling; but when Q£ — > 0, E g {x) tends 
to follow much closer to A{x) as expected. In addition, 
the minimum of E g {x) is always slightly higher than the 
minimum of A(x) by an amount that also diminishes 
upon Q£ — > 0. This behavior will be further clarified in 
the next section. 

The ratio E g /A vs. Uq/U or Q£ is plotted in FIG. 
ID The suppression of the gap strengthens when either 
the inhomogeneity becomes stronger (Uq/U is large) or 
its length scale L ~ l/Q becomes smaller, consistent 
the results in FIG. [3] The E g suppression relative to 
A, together with the fact the T^ F is largely deter- 
mined by strongest-coupling region, implies that the ratio 
2E g /Tl IF is generically reduced. The ratios 2E g /T* 1F 
and 2E 3 /Tkt are plotted in FIG. [5] for several repre- 
sentative cases. As expected, there is always a strong 
suppression of the ratio 2E g /T^ IF from 3.52; for a two- 
dimensional system, however, the ratios with Tkt are 
more subtle: for very small Q£ the ratio 2E 9 /Tkt might 
be enhanced due to the large deviation of Tkt from T//' IF 
(see also FIG. [2c)), while for larger value of Q£ the 
phase fluctuation region is narrow(see also FIG. [2)Ja)), 
and 2E 9 /Tkt is reduced from 3.52. 

For the purpose of comparison with the thin film ex- 
periments, a comment on the determination of T^ IF and 
Tkt is in order. Due to disorder and phase fluctuations, 
the resistive transition curve can be significantly broad- 
ened. T^ F can be estimated as the temperature at which 
the resistance drops to half of its normal state value, 
while Tkt can be defined as the temperature at which 
the resistance drops below the measurement threshold 
(see, for example, Ref. l28l ). Alternatively, one can 
extract T C MF from fitting the fluctuation resistance to 
Aslamazov-Larkin theory^ 2 -, and obtain Tkt from non- 
linear I-V characteristics or from fitting the resistance 
below T C MF to Halperin-Nelson formS (see, e.g., Refs. 
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1.8r 




00 

FIG. 3: The local order parameter A(x) and the local gap 
E g [x) (in units of A(Uq = 0) = 2u D e~ 1/UNF ) vs. spatial 
coordinate x £ [0,7r/Q]. Q£ = 0.3 and 0.1 in subfigure (a) 
and (b), respectively. UN F = 0.2, UqN f = 0.02. Here, U 
and Uq are the uniform and oscillating components of the 
coupling constant, respectively. N F is the density of states 
of the normal state; Q is the modulating wavevector of the 
inhomogeneous coupling constant; £ is the superconducting 
coherence length. The estimated numerical error of E g (x) is 
about 0.01. 

[HI!). Thus both T^ IF and T KT in principle can be 
measured from experiments, and can be used for com- 
parison with our theoretical results here. 



C. Additional inhomogeneities 

Apart from modulation of the coupling U, one may 
also be interested in a simultaneous modulation of other 
properties. For example, in the small Qt; limit, one may 
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00 

FIG. 4: The ratios of the spatially averaged gap E g to the 
spatially averag ed A (in units of 2Lu D e-^ /0NF ) (a) vs. Uq/U 
with Q£ = 0.3; (b) vs. Q£ with Uq/U = 0.1. UN F = 0.2 
in all cases. Here, U and Uq are the uniform and oscillating 
components of the coupling constant, respectively. Nf is the 
density of states of the normal state; Q is the modulating 
wavevector of the inhomogeneous coupling constant; £ is the 
superconducting coherence length. The estimated numerical 
error is about 0.01. 



expect the periodicity of U to be accompanied by a peri- 
odicity of the local density of states at the fermi level, or 
the mean free path. Another possible modulation, that 
of a periodic potential, is suggested in^, and in practice 
is equivalent to local modulation of U. Indeed, one may 
use an effective description of the self consistency equa- 
tion ([5]), taking Np — * Np + Nq cos(Qx) to lowest order 
in the amplitude Nq of the local DOS in the form: 

A(r) = N F U mod (r)TrTj2iUn (f)- (29) 

n 
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where U„ 



U 



NqU+N f Uq 
N f 



cos(Qx), and Nf is the 
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FIG. 5: The ratios of the spatially averaged gap E a to T c 
or T KT (a)vs. [/ Q /I7, Q£ = 0.3; (b) vs._U Q /U, Q£, = 0.1; 
(c) vs. Q£, (7q/C7 = 0.1. In all cases UN F = 0.2. Here, 
U and [/q are the uniform and oscillating components of the 
coupling constant, respectively. Nf is the density of states 
of the normal state; Q is the modulating wavevector of the 
inhomogeneous coupling constant; £ is the superconducting 
coherence length. The estimated numerical error is about 
0.02. 



spatially averaged DOS. Formally this is exactly the same 
as Eq. ([T]), and can be treated similarly, taking 



Uc 



N Q U + N F U Q 
N f 



(30) 



In practice, a local periodic potential may be imposed on 
the system externally by either acoustic means or an elec- 
tromagnetic field. Thus it might be interesting to check 
the change in Tc of a superconductor in the presence of 
an acoustic wave experimentally. 

Another possibility of interest is that along with U 
the electron mean-free path is modulated in the system. 
This can be naturally occurring if the periodicity in U is 
a consequence of spatial variation in the properties of the 
material used. Alternatively, one may obtain this case by 
a periodic doping of the superconductor. 

In this case we may describe the system effectively by 
modification of the Usadel equation © to: 

- -V- (DV0) = Acos6»-cj„sin6», (31) 

and taking the diffusion coefficient D to be spatially de- 
pendent. Choosing D = D + Dq cos{Qx) and repeating 
the treatment above, we find that Dq does not change the 
values of the Green's functions 9q,0i above (It however 
appears at higher orders of the equation) , and so doesn't 
change the results of this paper within this order. 



IV. SUPERCONDUCTOR-NORM AL-METAL 
(SN) SUPERLATTICE ANALOGY 

Some insight into the nature of the lowest-lying excita- 
tions for both large and small Q£ cases can be gained by 
considering a simplified system: superconductor-normal- 
metal-superconductor (SNS) junctions. First, consider a 
single SNS junction with length L = 2ir/Q, and A(x) = 
A, in the S, N part respectively. Andreev bound states 
will form in the normal metal, and the energy of these 
states can be obtained by solving Bogoliubov-de Gennes 
(BdG) equations for the clean case, or Usadel equations 
for the dirty case. In the limit L — * 0, the energy of 
the lowest-lying state is A, while in the opposite limit 
L ^> £, the (mini)gap is much smaller than A: in the 
clean case E g ~ vp/L ~ (Q£)A and in the dirty case the 
gap equals the Thouless energy D/L 2 ~ (QQ 2 ££L£££L. 
These states exponentially decay into the superconduc- 
tors for a distance ~ £. 

Based on a single SNS junction, one can build an SN 
superlattice with alternating superconductor and normal 
metal, each with length L = 2tt/Q, and A(x) = A, 
in the S, N part respectively. If L ^> £, Andreev bound 
states remain localized in the normal regions with the gap 
much smaller than A. On the other hand if L <C £, these 
states strongly mix with each other, and they form a 
tight-binding band. Therefore the gap, namely the lower 
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band edge, is lower than A, and in the limit Q£ — > oo it 
is precisely at A/2, the averaged A(x) (see the analytical 
calculation by Ref. [H) . The SN superlattice thus allows 
a qualitative understanding of the gap's behavior in the 
problem we addressed above: if Q£ ^> 1, all excitations 
are extended in space, with the uniform gap E g w A; if 
Q£ <C 1, the lowest- lying excitations are localized in the 
weakest coupling regions whose gap is close to the mini- 
mum of A(x). This analogy also elucidates the features 
in FIG. [3] given a point in space xq, E g (xo) is generally 
lower than A(xo), because the wave function of the low- 
lying excitations originating at a nearby region (within 
~ £) with smaller A(x) are exponentially suppressed at 
Xq, and when £ is smaller this effect is reduced; thus 
E g {x) follows closer to A(x) in the limit Q£ — > 0. Fi- 
nally, the difference between the minimum of E g (x) and 
the minimum of A(x) resembles the minigap in SN su- 
perlattice ~ vp/L or D/L 2 , which approaches zero as 
Q£ -> 0. 



SUMMARY AND DISCUSSION 



3.52 




FIG. 6: Schematic plot of the ratios 2E g /T* 1F and 2EJT K t 
vs. Q£. Here E g is the spatially averaged gap in local DOS; 
T C MF is the mean field T c ; Tkt is the Kosterlitz-Thouless 
transition temperature in 2d; Q is the modulating wavevec- 
tor of the inhomogeneous coupling constant; £ is the supercon- 
ducting coherence length. 1,2,3, and 4 are labels of different 
regimes described in the text. 



In this paper we investigated the properties of dirty 
BCS superconductors with a fluctuating pairing cou- 
pling constant U(x) = U + Uq cos(Qx). Particularly, 
we analyzed the change in the mean field T c , the zero- 
temperature order parameter A(x), and the energy gap 
in quasiparticle excitation E g (x) using the Usadel equa- 
tion for quasiclassical Green's functions. In addition, wc 
estimated the Kosterlitz-Thouless transition temperature 
Trt- Our analysis found four different regimes: 

(1) Q£, — > oo. In this case the mean field T c and the 
spatially averaged order parameter A are determined by 
the effective coupling constant U e ff > U [see Eq. (fT5|) ] . 
Moreover, since in this regime any quasiparticle wave- 
function is extended over the length scale L = l/Q, the 
local energy gap E g is uniform in space, and we found 
it to coincide with the spatially averaged A. The ratios 
2A/T C = 2E g /T c = 3.52 maintain their universal BCS 
value. 

(2) Q£ > 1. In this regime the physics is qualitatively the 
same as that of the previous case. The gap E g , however, 
is smaller than A by an amount that grows with decreas- 
ing Q£ or increasing Uq/U. Therefore 2E g /T c < 3.52 
(see FIG. HJ. 

(3) Q£ < 1. The system tends to divide into regions 
which behave according to the the local value of U(x). 
Thus the mean field T c is determined by the first for- 
mation of local superconductivity upon lowering temper- 
ature, and therefore T^ F is close to highest 'local T c \ 
In contrast, the global energy gap or the spatially aver- 
aged local gap is largely determined by the region with 
smallest U(x). Consequently, in this regime the ratio 
2E g /T^ IF is always suppressed from the universal BCS 
value, 3.52 (see FIG.Ek). Moreover, although the system 
is affected by phase fluctuations, in this regime Tjct is 
close to T* IF for small values of Uq (see FIG. Ha,). Thus 



2E g /T KT is also smaller than 3.52 (see FIG. [30 . 
(4) Q£ — > 0. As opposed to the previous regime, here 
phase fluctuations lead to a large suppression of Tkt rel- 
ative to T C MF (see FIG. Wp). Although 2E g /T™ F is still 
below 3.52, the ratios 2E 9 /Tkt is close to or larger than 
3.52 (see FIG. Et) L 

The value of 2E g /T* IF and 2E g /T KT vs. the entire 
range of Q£ is plotted schematically in FIG. [SJ with 
regimes 1-4 explicitly labeled in the graph. Schematic 
results of T C MF and Tkt vs. Q£ are summarized in FIG. 

m 

Finally, we discuss connections with thin film 
experiment s 48 ! 49 . A straightforward realization of in- 
homogeneous coupling is in disordered superconductor- 
normal-metal (SN) bilayer thin films. In a homogeneous 
bilayer SN with thickness smaller than the coherence 
length £, mean field analysis yields that T c and the en- 
ergy gap E g of the system are determined by the averaged 
coupling constan t 50 ! 51 ' 52 

where U is the pairing coupling in the superconducting 
layer, d is the thickness, N is the DOS at the Fermi 
energy, and the subscripts S and N denote the super- 
conductor and normal metal layers respectively. Thus 
the ratio 2E g ( T=0 * ) /T c is expected to remain at the BCS 
value 2it/C w 3.52 in a homogeneous SN bilayers. Never- 
theless, from (|32|) one observes that a spatially inhomoge- 
neous thickness cIs,n(x) (which is also consistent with the 
granular morphology of the sample^ 9 -) leads to a nonuni- 
form coupling U(x) even if the original coupling U is 
homogeneous. Therefore thickness variation generically 
leads to a superconductor with inhomogeneous pairing 
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FIG. 7: Schematic plot of the mean field transition tempera- 
ture T C MF and the Kosterlitz-Thouless temperature Tkt vs. 
Q£, where Q is the modulating wavevector of the inhomoge- 
neous coupling constant; £ is the superconducting coherence 
length; T c , ma x = i§-uj D e~ 1/NF(tJ+u Q ) is the maximum T C MF ; 
T c . a = ^-LUDe~ 1 / NFU is the mean field T c for a uniform cou- 
pling U. 1,2,3, and 4 are labels of different regimes described 
in the text. The qualitative feature of these results on T c are 
similar to those of Ref. [54| on clean superconductors. 



coupling. According to our results, a deviation of 2E g /T c 
from 3.52 is expected in such a system. 

Indeed our study was motivated by such observations. 
In Refs. I4al49l Long et al. report measurements of re- 
cently fabricated a series of Pb-Ag bilayer thin films, 
with thickness dpb = 4nm and d,A g increases from 6.7nm 
to 19.3nm. They observed a significant reduction of 
2E g /T^ IF from the expected value ~ 3.52, where E g 
is the spatially averaged gap extracted from tunneling 
measurement of the DOS, and T C MF is measured as the 
temperature at which R(T) drops to half of its nor- 
mal state value, and the resistive transition is sharp and 
well-defined. This suppression of 2E g /T^ IF is more pro- 
nounced in systems with thicker Ag thereby lower T C MF . 
In these samples with T^ F decreasing from 2.55K to 
0.72K with increasing d,Ag, the ratio 2 E a /T^ F decreases 
from - 3.6 to - 2.6 (see FIG 3(b) of Ref. El). 

These results can be qualitatively well understood by 
our study. The reduction of 2E g /T^ IF from 3.52, to- 
gether with the observed fact that the resistive transition 
is sharp and well-defined^ 8 , implies that the experimen- 
tal systems are in the regime (2) or (3) of our theoretical 
results summarized above (see FIG. [5]). In these regimes 
both 2E g /T^ F and 2E g /T KT are lower than 3.52, and 
the phase fluctuation is either absent or small enough 
to keep Tkt close to T^ F , explaining the sharp resis- 
tive transition. For samples with lower T c , U is smaller. 
Therefore, if we assume roughly the same amount of 
Uq for all samples, the effect of inhomogeneity will be 
stronger for samples with lower T c samples, and, conse- 
quently, the gap-to- T c ratio is even smaller for them. To 



3.5 




FIG. 8: The ratios of the spatially averaged gap E g to the 
mean field transition temperature T C MF or the Kosterlitz- 
Thouless transition temperature Tkt vs. UNf- UqNf = 
0.002, Q£ = 0.3. Here, U and Uq are the uniform and oscil- 
lating components of the coupling constant, respectively. Nf 
is the density of states of the normal state; Q is the modulat- 
ing wavevector of the inhomogeneous coupling constant; £ is 
the superconducting coherence length. Since T^ F monotoni- 
cally increases with U, this result resembles the experimental 
data of Ref. Es| (see FIG. 3(b) ibid.), which shows that the 
lower the measured T c of a thin-film bilayer is, the smaller the 
ratio 2_E 9 /T c .The estimated numerical error is about 0.02. 



make a rough comparison, we have calculated the gap-T c 
ratio vs. U for fixed Uq and plotted the results in FIG. [5] 
Although not claiming more than a qualitative explana- 
tion of the bilayer measurements, we note that our FIG. 
HI resembles FIG. 3(b) of Ref. EU 

An interesting venue for future research, which may 
extend to more 2d superconducting systems, is to con- 
sider a general fluctuation of the pairing interaction, not 
restricted to a particular wave number, but rather having 
a particular correlation length. In addition, aside from 
the low gap-T c ratio, Ref. [48] has also reported an unex- 
pected subgap density of states of quasiparticles in the 
same bilayer materials. Although our current model does 
not produce this behavior, one expects that it could be 
explained by including large spatial fluctuations of the 
pairing interaction (e.g. -fj- ~ 1), which strongly sup- 
press the gap, and the effect of mesoscopic fluctuations 
which tend to produce subgap states^. 
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APPENDIX A: CALCULATION OF A (T=0) IN 
THE LIMIT Q£ > 1 

Here we show some calculation details in deriving equa- 
tion (fTB|) . At T — the self-consistency equations are 



We take the limit x = ^ > 1 and a = (Q£) 2 > 1 
simultaneously, but their relative ratio might be cither 
large or small. Also using arctanh(z) = 1/2 ln(|l + z|/|l — 
z\), one can show that in this limit the above integral 
equals 



N F U 



du sin 8r 



N F U Q 



duiOi cos 9q 



Ai _ N F U 
~~2 ~ 2 



dw8\ cos 9q 



N f U q 



The evaluation of the integrals gives (define a 
and xq = a->£>/Ao): 



diu sin 6i 



DQ 2 /2 
A 



dio sin 9q = Aq arcsinh 



up 
Aq 



Aq In 



(2oj 



D 



\ A 



dioOi cos Oq 



Ai 
2a 



{—2 arctan(xo) + 2a arcsinh(a:o) 



- Va 2 -1 
+ arctanh ^ 

— 2 arctanh 



arctanh 



xoVa 2 ^ 1 + 1 



XoV a 2 — 1 — 1 



Ai 
2a 



Ai 
2 



2a ln(2xo) — a 
/ 2x^a 

I _a_ t £o _rr 

\ 2x n ^ 2a ^ 1 , 



1 / 2x a 

Z V 2:eo + 2a 1 . 



hi 



\x - a 



21n(2x ) 



In 



x + a 
2xqo, \ 



I Xq _ a 

I 2a 2x 



In 



= Ai In 1 + — = Ai In 1 



2u>p 
DQ 2 J UNp 



\xp - a\ 

Xq + a 
A! 



Ki, 



(AI) where K\ has exactly the same form as defined in (| 13[) 
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